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Abstract 

At the Large Hadron Collider, we prove the feasibility to detect pair production of 
the lightest CP-even Higgs boson /i of a Type II 2-IIiggs Doublet Model through 
the process ggC) Vhh (Higgs-strahlung, V = W^,Z), in presence of two h ^ bb 
decays. We also show that, through such production and decay channels, one has 
direct access to the following Higgs self-couplings, thus enabling one to distinguish 
between a standard and the Supersymmetric version of the above model: Xnhh ~ which 
constrains the form of the Higgs potential - as well as ^w^H^h ^ZAh ~ which are 
required by gauge invariance. Unfortunately, such claims cannot be extended to the 
Minimal Supersymmetric Standard Model, where the extraction of the same signals is 
impossible. 
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1 Introduction 



The ability to access Higgs self-couplings would enable one to constrain the form of the Higgs 
potential responsible for Electro- Weak Symmetry Breaking (EWSB) and, in particular, to 
assess whether such a potential corresponds to the one embodied in the Standard Model 
(SM) or indeed in extended versions of it. Lately, amongst the scenarios with an enlarged 
Higgs sector, much attention has been devoted to a (CP-conserving) Type II 2-Higgs Doublet 
Model (2HDM), possibly in presence of minimal Supersymmetry (SUSY) - the combination 
of the two yielding the so-called Minimal Supersymmetric Standard Model (MSSM) [1]. In 
a generic 2HDM and in the MSSM, of the initial eight degrees of freedom pertaining to the 
two complex Higgs doublets, only five survive as real particles upon EWSB, labelled as H, 
A (the first two are CP-even or 'scalars' (with Mh < Mh) whereas the third is CP-odd 
or 'pseudoscalar') and H^, as three degrees of freedom are absorbed into the definition of 
the longitudinal polarisation for the gauge bosons Z and W^, upon their mass generation 
after EWSB. This thus makes four additional Higgs particle states with respect to the SM, 
wherein only one complex Higgs doublet is allowed which, after EWSB, generates a lone 
CP-even (scalar) Higgs boson (henceforth also denoted by h). 

The detection of a h state, independently of the underlying model (amongst those men- 
tioned), is guaranteed at the Large Hadron Collider (LHC) [21 El S], over the theoretically 
allowed mass intervals. However, notice that while the Higgs boson mass Mh is a free param- 
eter in both the SM and a generic 2HDM, so that its upper limit (of about 700 GeV) is only 
dictated by unitarity requirements, in the MSSM the SUSY relations interlinking gauge and 
Higgs couplings imply an absolute upper limit on the mass of the lightest Higgs boson 
of about 140 GeV (see the second paper in [2] and references therein). Furthermore, there 
is a lower limit from LEP, of about 115 GeV, applicable directly to Mh in the SM but also 
adoptable in both a generic 2HDM and the MSSM over large areas of their parameter spaces. 
Finally, the Tevatron has put an exclusion limit on Mh in the SM in the range 162-166 GeV 

The CP-conserving 2HDM, with soft Z2 symmetry breaking, is identified uniquely at tree 
level by seven independent parameters (as will be made clear below): Mh, M^, M^, Mh±, 
tan/3 (the ratio between the vacuum expectation values of the two Higgs doublets), a (the 
mixing angle between the two CP-even neutral Higgs states) and A5 = (A5 — Aq)/2, one of 
the parameters of the Higgs potential as defined in [6] (see also [^). In the case of the MSSM, 
at tree level, only two parameters are needed, tan /3 and any of the Higgs masses (henceforth 
Ma)- (Recall that, in the SM, Mh is the only free parameter in the Higgs potential.) 

An intriguing situation that could arise at the LHC is then the following. That only a 
light Higgs state, with mass below 140 GeV, is found and it is impossible to establish whether 
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it belongs to the SM, a generic 2HDM or indeed the MSSM. In fact, this happens in the 
so-called 'decoupling region' of the extended scenarios, i.e., when Mh, Ma, Mh± 3> M/j (and 
for suitable choices of the additional free inputs in the generic 2HDM), quite independently 
of tan /3, where - for the same mass - the h couplings to ordinary matter in the SM are the 
same as those in both the 2HDM and MSSMj. Even in these conditions, however, it has 
been proved that one could possibly establish the presence of an extended Higgs sector by 
determining the size of the trilinear Higgs self-coupling Xhhh [8j, by looking at the Vector 
Boson Fusion (VBF) process, qq'^'^ — )■ qq^'^hh, and possibly also distinguish between a generic 
2HDM and the MSSM by isolating the Xnhh trilinear Higgs self-coupling in the same VBF 
channel . This can be achieved by looking at decays of the two light Higgs states produced 
into bottom-quark pairs. 

In general though, this distinction between a generic 2HDM and the MSSM is not always 
possible. Not even the discovery of also the A and states may say the conclusive word 
in this respect. In fact, despite there exist well established spectra among the four different 
masses in the MSSM (dictated indeed by SUSY for fixed, say, and tan/3), it may well be 
possible that the additional 2HDM parameters arrange themselves to produce an identical 
mass pattern. However, such a degeneracy between the two models would not typically 
persist if one were able to also measure other Higgs-gauge couplings, chiefly \iY±H^h and/or 
^ZAh, which are both proportional to cos(/3 — a) times gauge couplings. Even though the 
form of the Higgs-gauge coupling is specified by reasons of gauge invariance, the actual size 
can be different in the MSSM and in a generic 2HDM since in the latter both a and /3 
are independent parameters. Hence, the measurement of at least one of these Higgs-gauge 
couplings would constitute a confirmation, or otherwise, of the MSSM relations if one knew 
M/j (and possibly the heavier masses, M^, plus Mh±) and tan/3 but not a. 

These Higgs-gauge couplings are accessible via the double Higgs-strahlung process, that 
we intend to study hercl. Our approach is very conservative, as we make the assumption 
that only one parameter is known, M/j, measured for some value below 140 GeV, as may well 
happen at the LHC after only a h resonance is detected. We further imply that all SUSY 
states are much heavier than ordinary particles (with the possible exception of the lightest 
SUSY particle, which may well escape detection) and not yet detected by the time the h 
resonance is isolated and our advocated studies are pursued. In particular, in this paper, 
we show that, by resorting to double Higgs-strahlung, it is possible to extract both \w±H^h 



^In the MSSM, one also has to presume that the sparticle states are very heavy, so that the existence of 
a SUSY scenario is not obviously manifest through their detection. 

^In principle, they can also be accessed in VBF: see the last four graphs of Fig. 1 in [5^. In practise, we 
had verified in that paper that the contribution due to those diagrams was not accessible, the main reason 
being that they are not resonant in VBF, unlike the case of Higgs-strahlung. 
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and XzAh in a generic 2HDM whilst in the MSSM neither of them is accessible. 

Therefore, in conjunction with the measurement of Xuhh, accessible in VBF [9] and also 
in the process considered here (as we shall show), one will in principle be able to disentangle 
between the 2HDM and the MSSM, even if only a light Higgs resonance is found after earlier 
LHC data and no SUSY states are visible throughout. However, all this phenomenology 
requires that enough luminosity can be accumulated at the LHC, such that the Super-LHC 
(SLHC) option [TT] becomes needed to perform such kind of studies. 

We will illustrate how we have come to this conclusion, i.e., after investigating the pro- 
duction of lightest Higgs boson pairs in the Higgs-strahlung process, namely [111 [13] : 

g/) ^ Vhh, (1) 

with V = or Z and g*^'-* referring to any possible (anti)quark flavour. The relevant 
Feynman diagrams corresponding to this channel in both the 2HDM and MSSM considered 
here can be found in Fig. [T] (Herein, notice the three graphs at the bottom, where the Xhhh 
or Xuhh and Xw±H^h or XzAh couplings are located.) In our selection analysis, we will resort 
to the extraction of two h ^ bb resonances, in presence of both hadronic and leptonic V 
decays, i.e., we will study the following signatures: 

1. 'four 6-quark jets and two light-quark jets', emerging from ([1]) if W^, Z — )■ jet-jet; 

2. 'four 6-quark jets and a lepton pair', emerging from ([1]) if Z — t- 

3. 'four 6-quark jets, one lepton and missing energy', emerging from ([T]) if — )■ i^Ui. 

We will limit ourselves to the case i = e, fi and we will discuss the procedure to tag the 
6-jets. 

Our paper is organised as follows. In the next section, we outline the computational 
procedure. Sect. 3 presents our numerical results and discusses these. Sect. 4 contains our 
conclusions. 



2 Calculation 

We have assumed ^/s = 14 TeV for the LHC energy throughout. Our numerical results are 
obtained by setting the renormalisation and factorisation scales to 2Mh for the signal while for 
the QCD background we have used the average jet transverse momentum (p^ = J2iPTj/^)- 
Both Higgs processes and noise were estimated by using the Parton Distribution Function 
(PDF) set MRST99(COR01) [H]. All signal and background calculations were based on 

•^Similar conclusions also hold at linear colliders, as shown in Ref. (10) . 
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exact tree-level Matrix Elements (MEs) using either the ALPGEN program [15] or the 
HELAS subroutines [16]. Furthermore, VEGAS [17] or Metropolis [18j were used for the 
multi-dimensional integrations over the phase space. As for numerical input values of SM 
parameters, we adopted the ALPGEN defaults. 

Concerning the MSSM setup, the two independent tree-level parameters that we adopt are 
Ma and tan/3. Through higher orders, we have considered the so called 'Maximal Mixing' 
scenario (Xj = At — jj,/ tan/3 = a/6Msusy) [H]? wherein we have chosen for the relevant 
SUSY input parameters: /i = 200 GeV, Ai, = 0, with Msusy = 5 TeV, the latter - as already 
intimated - implying a sufficiently heavy scale for all sparticle masses, so that these are not 
accessible at the LHC and no significant interplay between the SUSY and Higgs sectors of 
the model can take placqj. Masses and couplings within the MSSM have been obtained by 
using the HDECAY program |2Q] . The tan /3 and Ma values adopted as representative of the 
MSSM parameter space were: tan /3 = 3 and 40 for Ma = 160, 200 and 240 GeV (yielding 
Mh,H = 108(178) GeV, Mh,H = 112(212) GeV and Mh,H = 114(249) GeV for the low tan/3 
value plus Mh,H = 129(160) GeV, M^^h = 129(200) GeV and Mh,H = 129(240) GeV for the 
high one, respectively). 

Before giving the details of the 2HDM setup we are using, let us recall the most general 
CP-conserving 2HDM scalar potential which is {Z2) symmetric under $1(2) ^(+)*&i(2) 
up to softly breaking dimension-2 terms (thereby allowing for loop-induced flavour changing 
neutral currents) [T], 

1 2 1 2 

V = ml^^l^i + m\^^\^2 - {mi2$i$2 + /^-c.} + -Ai ($l<l>i) + -A2 (^\^-^ + 

+ A3 ($l$i) (<l>|$2) + A4 ($1$ 2) + {^As (<I>1$2)' + /i.c.j . (2) 

In the following, the parameters mn, m22, Tfiyi-, Ai, A2, A3 and A4 are replaced by f , M^, M^, 
Mai Mh±, P and a (with v fixed). Hence, as intimated already, the CP-conserving 2HDM 
potential is parameterised by seven free parameters. Notice that from the scalar potential 
all the different Higgs couplings needed for our study can easily be obtained. (See [HJ [7] for 
a complete compilation of couplings in a general CP-conserving 2HDM.) In addition, since 
we are considering a so called Type II 2HDM the Yukawa couplings are also fixed by tan /3 
in the same way as in the MSSM. 

In our 2HDM, we will fix M^ and M^ to values similar to the ones found in the MSSM 
scenario we are considering, by adopting three different setups: 



^The only possible exception in this mass hierarchy would be the Lightest Supersymmetric Particle (LSP), 
as intimated already, whose mass may well be smaller than the lightest Higgs mass values that we will be 
considering. However, we have verified that invisible h decays (including the one into two LSPs) have 
negligible decay rates. 
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1. Mft = 115 GeV, Mh = 250 GeV, Ma = 220 GeV, Mh± = 300 GeV, 



2. Mh = 115 GeV, Mh = 300 GeV, Ma = 305 GeV, Mh± = 310 GeV, 

3. Mh = 115 GeV, Mh = 400 GeV, Ma = 405 GeV, Mh± = 410 GeV, 

where in the first scenario the masses have been chosen such that Mh > '2Mh, Ma > Mh+Mz 
and Mh± is above the hmit given by measurements of 6 — > S7 (for recent hmits see [2T| [22] 
and references therein). 

We always scan over the remaining parameters in the ranges 

-7r/2 < a < it/2, 
-471 < A5 < An, 
< tan/3 < 50. 

In order to accept a point from the scan we also check that the following conditions are 

fulfilled: the potential is bounded from below, the Aj's comply with the unitarity constraints 

of [23] and yield a contribution to |Ap| < 10"'^. In short, the unitarity constraints amount to 

putting limits on the eigenvalues of the S matrices for the scattering of various combinations 

of Higgs and EW gauge bosons. We have followed the normal procedure |I] of requiring the 

J = partial waves (ao) of the different scattering processes to fulfill |Re(ao)| < 1/2, which 

corresponds to applying the condition that the eigenvaluejf] Ay^_j_ of the scattering matrices 

or more precisely IQnS) fulfill \AyI±\ < Stt [2^. In other words we allow parameter space 

points all the way up to the tree-level unitarity constraint |Re(ao)| < 1/2. The spectrum of 

masses, couplings and decay rates in our 2IIDM is similar to the one exploited in Ref. pj^, 

obtained by using a modification of HDECAY [20j (consistent with a similar manipulation 

of the program used in Ref. [H]). For each accepted point in the scan the partial decay rates 

for the different Higgs bosons are then calculated using such HDECAY modification and also 

taking possible additional (i.e., 2HDM specific) partial widths of the H, A, and bosons 

into account. We have also cross-checked the results by using the 2HDMC program [12] 

to perform the same calculations. The results of the above discussed constraints on the 

parameter space are shown in Fig. [21 where the projections of the allowed parameters on 

three different planes. As-tan/?, sin(/3 — Q;)-tan/3 and sin(/3 — a)-A5, with a sample of 1000 

randomly chosen points are plotted. It is worth noting that the allowed region in the plane 

As-tan/? complies with the constraints found in Refs. [381 113], where the neutral Higgs- 

pair production and Higgs-strahlung at linear colliders within the general 2HDM have been 

investigated at the one-loop level. 

^Here, Z2 refers to the Z2 symmetry, Y is the hypercharge, and a is the total weak isospin. 

^For a variety of studies of Higgs pair production at the LHC, primarily through gluon-gluon fusion, see 

Refs. [lS]-[37]- Recent studies at linear and 77 colliders can be found in Refs. [55]) [32] and [ID]. 
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While the parameter dependence of the MSSM Higgs sector renders the computation 
of the tree-level MSSM cross-sections rather straightforward (as the latter depends on two 
parameters only, and tan the task becomes much more time-consuming in the context 
of the 2HDM. In order to calculate the cross-sections in this scenario, the latter are written 
as a combination of couplings and kinematic factors in the following way (n = 5): 



where all the explicit dependence on a, /3, Xuhh and Xhhh is contained in the couplings gi. 
gi = sin^(/3 - a), g2 = cos'^{(3 - a), g^ = cos(/3 - a)\Hhh, 9i = sin(/3 - a)\hhh, and g^ = 1, 
whereas the dependence on masses and other couplings is in the factors 



Note that the sum over subamplitudes Mj also contains all interference terms and that colour 
factors etc. are included properly. The cTjj are then calculated numerically for fixed masses. 
We can then get the cross-section in an arbitrary parameter space point by multiplying the 
kinematic factors with the appropriate couplings. However, there is a slight complication 
since the kinematic factor for the H hh contribution depends on the width Th if there 
is a s-channel resonance and the width in turn depends on the couplings. In this case the 
kinematic factor scales as 1/Th which is accounted for by assuming a fixed value for the 
width when the kinematic factor is calculated and then rescaling the result with the true 
width when calculating the contribution to the cross-section. The same type of procedure is 
also applied for the A — )■ Zh and — W^h resonant contributions. 

Notice that such a procedure, after the summation in eq. ([3]), amounts to the exact (i.e., 
complete and gauge invarianJll) calculation of the cross section for process ([1]). 

3 Results 

3.1 Inclusive Signal Results 

In this section, after a preliminary analysis of the Higgs mass and coupling spectra in the 
MSSM and a general Type II 2HDM, we will start our numerical analysis by investigating 
the model parameter dependence of the Higgs pair production process in ([1]) at fully inclusive 
level, assuming on-shell production of gauge and Higgs bosons, followed by the decay of the 
former in all possible final states (of quarks, leptons and neutrinos, as appropriate) and of 

^ Apart from subleading effects of 0{T/M), where M and F are the mass and width, respectively, of a 
resonant particle. 




(3) 




(4) 
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the latter into two bb pairs, with the integration over the phase space being performed with 
no kinematical restrictions. This will be followed by an analysis of the production and decay 
processes pertaining to the Higgs signals at fully differential level, in presence of detector 
acceptance cuts and kinematical selection constraints. Finally, we will compare the yield 
of the signals to that of the corresponding backgrounds and perform a dedicate signal-to- 
background study including an optimisation of the cuts in order to enhance the overall signal 
significance. We will treat the MSSM and 2HDM in two separate subsections. 



3.1.1 MSSM 

As representative of the low and high tan /3 regime, we will use in the remainder the values 
of 3 and 40. We have instead treated Ma as a continuous parameter, varying between 100 



and 700 GeV or sc^. Fig. [3] presents the fully inclusive MSSM cross-section for the process 
of interest, as defined in ([T]), followed by two h bb decays and of those of the and 
Z bosons into anything. In the same figure, we also show the rates corresponding to the 
case in which we have removed from the list of diagrams the last two in Fig. in order 
to individuate the origin of possible dynamic structures. In fact, the difference between the 
red and black lines illustrate the effects of the diagrams involving the Xw^H'fh and XzAh 
couplings, which can amount to a resonant contribution (at low tan /3) or a constructive 
interference (at large tan/3). 

The displayed rates however correspond to the ideal situation in which all final state 
objects (quarks and gauge bosons) are detected with unit efficiency and the detector coverage 
extends to their entire phase space, so that they only serve as a guidance in rating the 
phenomenological relevance of the processes discussed. 

A more realistic analysis is in order, which we have performed as follows. In all channels, 
the four 6-jets emerging from the decay of the hh pair are accepted according to the following 
criteria: 

p^T > 30 GeV, Ir^^l < 2.5, ARbb > 0.7. (5) 

Their ta ggin g efficiency is taken as eb = 50% for each b satisfying these requirements, e;, = 
otherwisqlj. We assume no 6-jet charge determination. In addition, to enforce the recon- 
struction of the two Higgs bosons, we require all such 6's in the event to be tagged and that 
at least one out of the three possible double pairings of 6-jets satisfies the following mass 



Values of Ma below 90 GeV or so are actually excluded by LEP for the lower tan /? value: see 01] . 
^Notice in fact that the two sets of diagrams peak in different phase space regions, so that, modulo the 
aforementioned 0{T/M) effects, they can safely be separated. 

-'^°Here and in the remainder, the label b refers to jets that are 6-tagged while j to any jet (even those 
originating from 6-quarks) which is not. As probability for a gluon or light-quark jet to fake a b jet we take 
R= 1/100. 
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preselection: 

imb,,b, - Mhf + (mb3,fe, - Mhf < 2 a^, (6) 

where a^n = 0.12 whereas 6162 and 6364 refer, hereafter, to the two pairs of 6-jets best 
reconstructing the Higgs masses with hi being the 6-jet with highest pt- (As intimated, we 
will be working under the assumption that a light Higgs state h is found already.) 

As for tagging the gauge bosons produced via the Higgs-strahlung process ([T]), there are 
two possibilities, depending on whether the gauge boson decays are 'leptonic' or 'hadronic'. 
In the former case, the additional cuts are 

p'j^'' > 20 GeV, It]'^"^'''] < 2.5, ARf,e,bmiss > 0.4, |M,+,- - Mz\ < 5 GeV (7) 

(with i = e, fi and wherein the neutrino transverse momentum is equated to the missing one, 
p™"^*^, while the longitudinal one is reconstructed by assuming M^j^^ = M\y± and choosing 
the solution with smallest magnitude). In the latter case, we use 

p^T>20GeV, \r]^<2.5, ARbjjj > 0.7, 70 GeV < A% < 100 GeV, (8) 

where Mjj is the invariant mass of the two non-6-jets in the event. 

In our investigation of hadronic final states, we will assume that 6-quark jets are distin- 
guishable from light-quark and gluon ones and, as already mentioned, neglect considering 
6-jet charge determination. In the case of leptonic ones, we will limit ourselves to the case 
of electrons and muons (collectively denoted by the symbol i). Finite calorimeter resolu- 
tion has been emulated through a Gaussian smearing in transverse momentum, p^, with 
ia{pT)/pTy = (0.60/y^)2 + (0.04)2 jg^g {(y{pT)/pT? = (0.12/^)2 + (0.01)^ 

for leptons. The corresponding missing transverse momentum, p™*^*^, was reconstructed from 
the vector sum of the visible momenta after resolution smearing. Also, in our parton level 
analysis, we have identified jets with the partons from which they originate and applied all 
cuts directly to the latter, since parton shower and hadronisation effects were not included 
in our study. 

Tab. [T] shows the rates of the signal after the implementation of the acceptance and 
preselection cuts in eqs. (IS])-® (hereafter, also referred to as 'primary cuts'). Notice that 
6-tagging efficiencies are not taken into account in this table. In fact, for the '46-jet' tagging 
option that we are advocating, one should multiply the numbers in Tab. [Uby ef, that is, 
1/16. 

The conclusion here is rather straightforward: given the smallness of the inclusive signal 
rates, the scope of extracting hh — )■ 46 signals in the MSSM from the double Higgs-strahlung 
process ([T]) is basically nil. Therefore, in our forthcoming signal-to-background analysis, we 
will not pursue the MSSM case. 
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Ma (GeV) 


(T(gg(') Vhh) (fb) 


160 


0.028(0.034) 


200 


0.019(0.037) 


240 


0.022(0.053) 


tan (3 = 3 


160 


0.013(0.031) 


200 


0.011(0.024) 


240 


0.0087(0.020) 


tan /3 = 40 



Table 1: Cross-sections for Higgs-strahlung after Higgs/gauge boson decays (relevant Branch- 
ing Ratios (BRs) are all included) and the acceptance and preselection cuts defined in (I5])-(|8]), 
for two choices of tan /3 and a selection of values, assuming the MSSM in Maximal Mixing 
configuration. Here, the numbers outside (inside) brackets refer to the 'leptonic'('hadronic') 
decay channel of the gauge boson as described in the text. No 6-tagging efficiencies are 
included here. 

3.1.2 2HDM 

As already alluded to earlier the parameter space of the general CP-conserving Type II 
2HDM we are considering here is quite large as it depends on seven independent parameters. 
In order to get a feeling for the dependence of the signal cross-section of process ([1]) upon 
them, we present in Fig. H] the results of our selected mass scenarios where we scan the 
allowed parameter space over 1000 randomly chosen points. 

Comparing with the cross-sections obtained in the MSSM we notice that in the 2HDM, 
see Fig. HI they can be up to two orders of magnitude larger. The main differences are due 
to the fact that: 

• the trilinear Higgs couplings are not restricted in size by SUSY relations; 

• the different parameters can vary independently of each other. 

At the same time the kinematic factors will be the same for a given set of masses and widths 
of the different Higgs bosons. Therefore, in those cases many features of the signal, such as 
differential distributions, will be similar to the MSSM even though the normalisation can 
be completely different. We also see from Fig. Hlthat there is quite a strong dependence on 
the masses of the Higgs bosons, with the maximal cross-section essentially decreasing by one 
order of magnitude as the heavy Higgs boson masses are increased from the scenario with 
lowest masses considered to the one with the highest masses. 
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As it is clear from Fig. |5l the largest inclusive cross-sections are obtained from the 
resonant channel H — )■ hh, whereas the contributions from the A — )■ Zh and — )■ W^h 
channels at the most reach about 10-20 fb. It is also clear from the figures that, even 
though there are slight differences in the parametric dependencies of the three different 
resonant contributions, the overall picture is very similar in all cases. 

Finally, we have also verified, though not shown, that the most promising signals are 
the ones where the decays leptonically, which is sensitive to both the H"^ — )■ W^h 
and H hh resonant enhancements, and possibly also the leptonic decay of the Z, which 
is sensitive to possible A — t- Zh and H hh resonances, whereas it is evident that the 
hadronic decay channels will not be relevant phenomenologically since they are only a factor 
~ 2-5 larger than the leptonic ones, yet burdened by overwhelming QCD backgroundJ"]. 

3.2 Signal-to-Background Differential Analysis 

The backgrounds that we have considered are as follows. In the case of process ([T]) for 
V = W^, they are: 

1. Wbbbb, irreducible (i.e., it has exactly the same particle content as the signal); 

2. tt, with two mistags (i.e., when two light- quark jets are erroneously classified as 6-jets); 

3. Wbbjj, with two mistags; 

4. Wjjjj, with four mistags; 

5. Whbb — Wbbbb, irreducible and co-resonant (i.e., it resonates at the h mass, albeit 
only once, unlike the signal, which does so twice); 

6. Whjj — )■ Wbbjj, CO- resonant with two mistags. 
In the case of process ([T]) for V = Z, they are: 

1. Zbbbb, irreducible; 

2. Zbbjj, with two mistags; 

3. Zjjjj, with four mistags; 

4. Zhbb —7- Zbbbb, irreducible and co-resonant; 

5. Zhjj — !■ Zbbjj, co-resonant with two mistags. 

^^Hence, we will not consider hadronic signatures any further in this study. 
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In all cases both QCD and EW channels were considered, as appropriate. The former by 
exploiting ALPGEN and the latter the HELAS subroutines (including all 2HDM background 
channels to the signal but not the interferences). We have however verified that, already 
after the acceptance and preselection cuts, the last two channels (5. and 6. for V = 
plus 4. and 5. for V = Z) are negligible in both cases, so that we have ignored them in the 
reminder of the analysis. In fact, also the backgrounds 4. for V = and 3. for V = Z are 
very small so that they could be neglected in first approximation. We have however included 
them in our simulation. The backgrounds 3. for V = and 2. ioi V = Z are similar both 
in magnitude and shape to the corresponding irreducible ones (i.e., 1. for both V = and 
V = Z) after taking into account the appropriate weights, including the 6-tagging/rejection 
factors discussed previously, although for clarity we do not plot them. Further, notice that 
background 2. for V = is significant before the final selection cuts, so that we have 
included it in our simulations. 

We start by investigating the differential structure of process ([T]) in the attempt to extract 
the contributions that are most sensitive to the Xuhh^ ^w^H^h = cos(/3 — a), and \zAh = 
cos(/3 — a) couplings, respectively. (Unfortunately, we can anticipate already that our studies 
revealed impossible to extract the h hh component of process ([T]).) To this end, we 
isolate in ([3]) the corresponding "resonant" contributions, more specifically a^z for Xuhh and 
(T22 for \w±H^h and A^^J^- Note that this corresponds to the approximation that the 
respective final states are WH or ZH and H^h or Ah, with the consequent decays H — )■ hh, 
— W^h and A — > Zh, for the two different processes at hand in ([1]). We have then 
identified, for the three different mass scenarios 1.-3. introduced in Sect. 2, those differential 
spectra that show significant differences from the yield of all the backgrounds above. 

In Figs. EHHl the resulting distribution^ after the primary cuts in eqs. ([S])-([H]) for the 
leptonic and Z decays are given in the 2HDM with heavy Higgs masses in the range 
from 200 GeV to ~ 400 GeV, limited to the contributions from the resonances, H ^ hh 
(with the H produced in association with a W^) in Fig. El H"^ — )■ hW^ in Fig. [3, A — > Zh 
in Fig. [H] plus H hh (with the H produced in association with a Z) in Fig. O Note 
that in Fig. [TJ as well as in the following Figs. [10] and [12] and Tab. H] only the two mass 
scenarios with Mh± = 310, 410 GeV are shown since the results for the Mh± = 300 GeV 
scenario is very similar to the one for Mh± = 310 GeV. The plotted observables are as 
follows: M^'° - the minimal mass of any pair of 6-jets, M^^,™™ - the next-to-minimal mass of 
any pair of 6-jets, p!^^ - the transverse momentum of the reconstructed h containing the 6-jet 



^^We have verified, through BRS tests, that, ahhough such contributions are not separately gauge invari- 
ant, gauge violating effects always occur at the permille level over the kinematic regions we are interested in 
(again, modulo the usual 0{T/AI) effects). 

^•^Notice that the normalisation in these figures is to unit area. 
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with highest transverse momentum, ) - the transverse momentum of the reconstructed 

W{Z) boson. In general terms, in all cases, the larger the Higgs masses involved in the decay 
the more effective the use of the described observables for the isolation of the signals. 

The shape of the plotted distributions, after accounting for correlations amongst the 
kinematic variables, leads to the following choice of "optimised cuts" in order to enhance 
the signals of interest compared to the backgrounds trying to isolate the couplings Xuhh, 
Xw^H'fh and XzAh 

• Xnhh from qq HW with H hh 

M^"' > 100 GeV, M"'"'" > 100 GeV, > GeV, > 100 GeV (9) 

• Avi/±_ff=F/i from qq — )■ H^h with — )■ hW^ 

M^'"" > 100 GeV, M^7^'° > 100 GeV, p^' > 100 GeV, > GeV (10) 

• XzAh from qq — > Ah with A — hZ 

M™" > 70 GeV, M°™" > 100 GeV, p!}.' > 100 GeV, p| > GeV (11) 

• Xnhh from qq — )■ HZ with H ^ hh 

Mjf" > 40 GeV, MX"' > 70 GeV, p^' > GeV, p| > 100 GeV (12) 

The resulting cross-sections for the relevant backgrounds after applying these optimised 
cuts are given in Tab. [2l After taking into account the 6-tagging {e^ = 0.5) and mis-tagging 
{R = 0.01) factors, this means that only the Vbbbb and Vbbjj backgrounds are important, 
although in our numerical analysis we have kept all the discussed backgrounds. 

In the next step of our analysis, we have scanned again the 2HDM parameter space in 
our three mass scenarios to search for points where, in presence of the above optimised cuts, 
one can establish a resonant signal at the 5a level. Since the number of expected background 
events {B) is very small we have used Poisson statistics in way similar to the one suggested 
in [l5] and calculated the number of signal events 5* that gives a probability for detecting 
S* + i? or more events, given B expected ones, to be less than 2.8 ■ 10~^ after / Cdt = 3000 
fb^^ of integrated luminosity, thus corresponding to a 5a significance. The required signal 
cross-sections obtained in this way are given in Tab. [31 

The results of the scans are given in Figs. [10] and [Til From looking at the latter, there is 
obviously scope in finding a 5a signal in almost all cases, for both the charged and neutral 
current version of process ([Ijll^, as shown in Figs. [TDland[TTl respectively. 



i-^The only notable exception is the hW^ resonance for Mh = 250 GeV, Mh = 115 GeV, Ma = 220 

GeV and M^± = 300 GeV. This is due to the cut M^^"' > 100 GeV, which is necessary against the tt 
background but at the same time kills the signal in this case. 
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cut-comb 


aipp — )■ Wbbbb) 


a(pp — i- Wbbjj) 


(j(pp — i- Wjjjj) 


a(pp -^tt^ Wbbjj) 




(fb) 


(fb) 


(fb) 


(fb) 


© 


9.2 ■ 10-4 


1.4 


82 


0.32 


(HOD 


1.8 • 10-3 


2.9 


1.6 ■ 10^ 


0.59 


cut-comb 


(7(pp — )■ Zbbbb) 


aipp — )■ Zbbjj) 


aipp — )■ Zjjjj) 

\rr J J J J / 






(fb) 


(fb) 


(fb) 




m 


3.9-10-4 


1.1 


19 




m 


4.0 ■ 10-4 


1.2 


20 





Table 2: Cross-sections for the most important background processes after the optimised 
cuts defined in fl9|)- f|T2l) . respectively, have been applied in addition to the acceptance and 
preselection cuts defined in (|5])-(|H]). Note that 6-tagging (e?, = 0.5) and mis-tagging [R = 
0.01) factors have not been applied. 

As it is clear from the plots the cross-sections with a — )■ W^h or a A — )■ Zh resonance 
are typically dominating the total cross-sections after the optimised cuts whereas the H hh 
resonance is typically not dominant even after applying them. The reason for this is that the 
width of the A and H"^ are typically small in the 2HDM scenarios we are considering and 
thus these resonances always contribute significantly to the total cross-section. In contrast, 
the H ^ hh is typically wide and it cannot be small at the same time as the A or are 
wide. 

Finally, Fig. [12] shows the phenomenal level of purity of the resonant signal samples that 
can be achieved, after the optimised cuts are implemented, for the "best case scenarios" in 
the 2HDM as given in Tab. IH Note that here the signals have been calculated from the sum 
of all diagrams in and thus all interferences are accounted for whereas the scenarios have 
been chosen to give the highest resonant cross-sections (the short-dashed lines in Figs. [10] 
and [TTj) . The high purity is exemplified by the invariant masses that would enable the 
reconstruction of the decaying Higgs state, i.e., the four 6-jet invariant mass (for the case 
of H ^ hh) and the invariant mass of the Higgs h not containing the highest transverse 
momentum 6-jet with the lepton-lepton (for the case A — )■ hZ) and lepton-neutrino (for the 
case if^ — )■ hW^) pair. (Here, the ti background is not shown, as it has become totally 
negligible.) As a final note we have investigated how far the different "best case scenarios" in 
Table [1] are from the strong coupling limit. By applying the coupled renormalisation group 
evolution equations for the A-couplings in Eq. (|2|) as well as the t, 6, r Yukawa couplings and 
the gauge couplings jl6] , we can determine how far a given point in parameter space can be 
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cut-comb 
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S 


a{qq Whh Ivibhbh) (fb) 
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0.30 


6.1 


0.033 


(HDD 


0.59 


7.3 


0.041 


cut-comb 


B 


S 


a{qq^ Zhh iibbbb) (fb) 


m 


0.16 


5.1 


0.027 


m 


0.16 


5.1 


0.027 



Table 3: The number of background {B) events expected after / Cdt = 3000 fb~^ of inte- 
grated luminosity when applying the optimised cuts defined in fP)- f[T2]) (in addition to the 
acceptance and preselection cuts) and the number of signal (S) events that give a probability 
less than 2.8 • 10^^ to observe S + B events when expecting B events using Poisson statis- 
tics. In addition the corresponding signal cross-sections are given (without any 6-tagging 
efficiency) . 

evolved in renormalisation scale before it reaches the tree-level unitarity limit This in turn 
gives an indication at what scale one can expect to see additional states and/or interactions 
(a.k.a. UV-completion). Applying this procedure we find that for most of the "best case 
scenarios" , the tree-level unitarity limit is reached at a scale which is a factor 3-4 larger than 
the input scale. The only two exceptions are the two scenarios with sin(/3 — a) = —0.01 in 
Table H] where this limit is reached only after evolving to a scale which is a factor 15 — 30 
times larger than the input scale. 



4 Conclusions 

In summary, we have shown the potential of a high luminosity LHC in extracting (resonant) 
decay signals of heavy Higgs states, both charged [H"^) and neutral {H and A), in a generic 
CP-conserving Type II 2HDM, using leptonic decays of a gauge boson and 6-quarks ones of 
the lightest Higgs state of the model (h), the latter produced in the Higgs-strahlung process 
gg*-'-* — )■ Vhh, with V = or Z and g*-'-* referring to any possible (anti)quark fiavour. The 
expanse of parameter space of the enlarged Higgs model that can be covered is sizable, not 
fine-tuned and yields detectable and extremely pure signals. This opens the prospect of 
accessing fundamental triple (Higgs and gauge) couplings of the model, notably Xuhh (but 
not Xhhh) plus Xw^H^h and XzAh, which would enable one to distinguish between this Higgs 
scenario and alternative ones (such as the MSSM) even in presence of degenerate mass and 
coupling spectra between the two. The possibility of extracting the latter and the precision 

-•^^We would like thank Oscar Stal for providing us with a code for solving the RGEs. 
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Mh = 300 GeV, M 


^ = 305 GeV, Mh± = 310 GeV 




{qq HW hhW) 


-0.26 


-2.9 


0.9 


-988 


387 


^ {qq H^h hhW) 


-0.01 


-1.3 
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-335 


-873 


(HH) {qq ^ Ah^ hhZ) 


-0.01 


-1.4 


2.9 


-397 


161 


(US]) {qq ^ HZ ^ hhZ) 


0.09 


-1.3 


24 


-992 


-293 


Mh = 400 GeV, Ma = 405 GeV, Mh± = 410 GeV 


dH]) {qq HW hhW) 


0.73 


-0.7 


5.2 


867 


-1328 


(HDD {qq H^h hhW) 


0.0034 


-2.5 


7.5 


-671 


68 


(dH) {qq Ah-y hhZ) 


-0.0012 


-2.5 


6.4 


-651 


58 


(jni) {qq ^ HZ ^ hhZ) 


0.73 


-0.7 


5.2 


867 


-1328 



Table 4: The respective best case scenarios used for illustration in Fig. [121 after the im- 
plementation of the optimised cuts defined in (l9i)- (fT2l) . Note that there is no entry for the 
H^ resonant scenario for Mh± = 300 GeV since the results are very similar to the ones for 
Mh± = 310 GeV. 

at which this can be achieved, depending on the parameter space regions considered, ought 
to be assessed in a full detector environment, in presence of all necessary parton shower and 
hadronisation effects. However, the level of sophistication of our study and the encouraging 
results obtained hint at the feasibility of these analyses. To confirm this, we make available 
upon request all our codes. 
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Diagrams by FEYMANDRAW 




Figure 1: Feynman diagrams for qiq2 — > V^h^h^. 
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2HDM: mi,=115 GeV, mjj=250 GeV, m^=220 GeV, mjj+=300 GeV 
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2HDM: m|^=115 GeV, mjj=400 GeV, m^=405 GeV, mjj^=410 GeV 
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Figure 2: The projection of the allowed points in the 2HDM parameter space on three 
different planes, As-tan^, sin(^ — Q;)-tan^ and sin(^ — a)-X5, in the three different mass 
scenarios: (top) Mh = 250 GeV, = 115 GeV, Ma = 220 GeV and M^± = 300 GeV, 
(middle) Mh = 300 GeV, Mh = 115 GeV, Ma = 305 GeV and Mh± = 310 GeV, (bottom) 
Mh = 400 GeV, Mh = 115 GeV, Ma = 405 GeV and Mh± = 410 GeV. The projections 
have been obtained with a sample of 1000 random points in the allowed 2HDM parameter 
space. 
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pp ^ V h h ^ jj/11 4b (E,^ = 14 TeV) 
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Figure 3: The inclusive cross- sect ions (as defined in the text) for Higgs-strahlung in ([T]), 
followed by hh — )■ bbbb and V — )■ 'anything' decays, as a function of the CP-odd Higgs boson 
mass, Ma, for two choices of tan/3, assuming the MSSM in Maximal Mixing configuration. 
The red lines correspond to the case in which we have removed from the computation the 
last two Feynman diagrams in the previous figure. 
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pp Vhh Xbbbb in 2HDM, mj,=l 15 GeV, mj^=250 GeV, m^=220 GeV, mjj^=300 GeV 
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pp Vhh Xbbbb in 2HDM, mj^=l 15 GeV, mj^=400 GeV, m^=405 GeV, m^^=4l0 GeV 
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Figure 4: The inclusive cross-section in the 2HDM (as defined in the text) for Higgs-strahlung 
in ([T]), followed by hh — )■ bbbb and V — )■ 'anything' decays in the three different mass scenarios: 
(top) Mh = 250 GeV, M/, = 115 GeV, Ma = 220 GeV and Mh± = 300 GeV, (middle) 
Mh = 300 GeV, M;, = 115 GeV, Ma = 305 GeV and Mh± = 310 GeV, (bottom) Mh = 400 
GeV, Mh = 115 GeV, Ma = 405 GeV and Mh± = 410 GeV, as a function of the remaining 
three parameters, for 1000 random points in the allowed 2HDM parameter space. 
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pp ^ VH ^ Vhh Xbbbb in 2HDM, m^=l 15 GeV, mj^=300 GeV, m^=305 GeV, mjj^=310 GeV 
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Figure 5: The contributions to the inclusive cross-section in the 2HDM (as defined in the 
text) for Higgs-strahlung in ([T]), followed by hh — bbbb and V — )■ 'anything' decays from the: 
(top) H — i- hh, (middle) A Zh, and (bottom) — )■ W^h, channels respectively in the 
mass scenario: Mh = 300 GeV, Mh = 115 GeV, Ma = 305 GeV and Mh± = 310 GeV, as 
a function of the remaining three parameters, for 1000 random points in the allowed 2HDM 
parameter space. 
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Figure 6: Differential distributions (normalised to unit area) of the H ^ hh contribution to 
Higgs-strahlung (solid lines) after the primary cuts of eqs. ([5])-(l8]) in the 2HDM (as defined 
in the text) followed by hh — )■ hhhb and — )■ iv^ decays in the different mass scenarios 
indicated. Background spectra from the Wbbbb and tt processes are also given (crosses). 
The distributions are as follows: the minimum bb invariant mass (top left); the next-to- 
minimum bb invariant mass (top right); the transverse momentum of the reconstructed Higgs 
boson containing the 6-jet with largest transverse momentum (bottom left); the transverse 
momentum of the reconstructed (bottom right). 
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Figure 7: Differential distributions (normalised to unit area) of the — )■ hW^ contribution 
to Higgs-strahlung (solid lines) after the primary cuts of eqs. (j5])-([8j) in the 2HDM (as defined 
in the text) followed by hh bbbb and — )■ iui decays in the different mass scenarios 
indicated. Background spectra from the Wbbbb and ti processes are also given (crosses). 
The distributions are as follows: the minimum bb invariant mass (top left); the next-to- 
minimum bb invariant mass (top right); the transverse momentum of the reconstructed Higgs 
boson containing the 6-jet with largest transverse momentum (bottom left); the transverse 
momentum of the reconstructed (bottom right). 
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Figure 8: Differential distributions (normalised to unit area) of the A hZ contribution to 
Higgs-stralilung (solid lines) after the primary cuts of eqs. ([5])-(|8]) in the 2HDM (as defined in 
the text) followed by hh —¥ bbbb and Z £i decays in the different mass scenarios indicated. 
Background spectra from the dominant irreducible background Zbbbb are also given (crosses). 
The distributions are as follows: the minimum bb invariant mass (top left); the next-to- 
minimum bb invariant mass (top right); the transverse momentum of the reconstructed Higgs 
boson containing the 6-jet with largest transverse momentum (bottom left); the transverse 
momentum of the reconstructed Z (bottom right). 
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Figure 9: Differential distributions (normalised to unit area) of the H ^ hh contribution to 
Higgs-strahlung (solid lines) after the primary cuts of eqs. ([5])-(|8]) in the 2HDM (as defined in 
the text) followed by hh — )■ hbhh and Z It decays in the different mass scenarios indicated. 
Background spectra from the dominant irreducible background Zbbbb are also given (crosses). 
The distributions are as follows: the minimum bb invariant mass (top left); the next-to- 
minimum bb invariant mass (top right); the transverse momentum of the reconstructed Higgs 
boson containing the 6-jet with largest transverse momentum (bottom left); the transverse 
momentum of the reconstructed Z (bottom right). 
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Figure 10: Distribution of resulting cross-sections in the 2HDM with leptonic decays 
when scanning over 10000 parameter space points for the different mass scenarios indicated. 
The two columns show the resulting cross-sections after applying optimal cuts (as defined 
in the text) to select the — )■ hW^ (left) or H ^ hh (right) channels respectively. The 
different lines show: the inclusive cross-section (solid), the cross-section corresponding to the 
— )■ hW^ (left) OT H ^ hh (right) resonances (dashed blue), and the same cross-section 
when requiring that the contributions from the resonance is at least 50% (red long dash). 
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Figure 11: Distribution of resulting cross-sections in the 2HDM with leptonic Z decays 
when scanning over 10000 parameter space points for the different mass scenarios indicated. 
The two columns show the resulting cross-sections after applying optimal cuts (as defined 
in the text) to select the A — )■ hZ (left) or H ^ hh (right) channels respectively. The 
different lines show: the inclusive cross-section (solid), the cross-section corresponding to 
the A — )■ hZ (left) ot H ^ hh (right) resonances (dashed blue), and the same cross-section 
when requiring that the contributions from the resonance is at least 50% (red long dash). 



> 

ID 

o 



^ 0.002 - 
b 



dashes: M4 = 230 GeV 
dots: Mj = 305 GeV 
dotdash: Mj = 405 GeV 
w Zbbbb 



> 



S 0.0010 
b 



k-<f yj'flit'y«^)cl»f")<-j < "' - ^"-Nr') j 'y-K')<-»v - L -!i « It' I 

500 600 

M(b3b411) [GeV] 



000 ^ — — 'tf7c^w^w»<-><-j < ' - a'y) g 'yw^w -X ' V - fa -i k )t ■ «! w 0000 

200 300 400 500 600 700 



dashes: = 250 GeV 
dots: Mh = 300 GeV 
dotdash: M„ = 400 GeV 
" Zbbbb 



200 300 



> 



^ 0.004 - 



T 



T 



X) 



b 



dots: M„+ = 310 GeV 
dotdash: Mj,+ = 410 GeV 
K ¥bbbb 



> 
CD 
O 



S 0.0010 
b 



[ h * n ri ) 



400 500 600 

M(4b) (GeV) 



dashes: Me = 250 GeV 
dots: M„ = 300 GeV 
dotdash: M„ = 400 GeV 
" Wbbbb 



M(b3b'lli/) [GeV] 



M('lb) (GeV) 



Figure 12: Differential distributions (normalised to the cross-sections in the respective "best 
case scenarios" ) after the optimised cuts for the four different signals: pp — > Ah with A hZ 
(top left); pp — )■ ZH with H hh (top right); pp H^h with — )■ hW^ (bottom left); 
pp — )■ WH with H ^ hh (bottom right). The distributions are as follows: the invariant mass 
of the reconstructed A boson using the second reconstructed h (top left); the invariant mass 
of the reconstructed H boson (top right); the invariant mass of the reconstructed H"^ boson 
using the second reconstructed h (bottom left); the invariant mass of the reconstructed H 
boson (bottom right). 



